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Abstract 


We present the results of a Monte Carlo technique to calculate the absolute magnitudes {H) and slope 
parameters (G) of ~ 240,000 asteroids observed by the Pan-STARRS 1 telescope during the first 15 months 
of its 3-year all-sky survey mission. The system’s exquisite photometry with photometric errors < 0.04 mags, 
and well-defined filter and photometric system, allowed us to derive accurate H and G even with a limited 
number of observations and restricted range in phase angles. Our Monte Carlo method simulates each 
asteroid’s rotation period, amplitude and color to derive the most-likely H and G, but its major advantage 
is in estimating realistic statistical-Hsystematic uncertainties and errors on each parameter. The method was 
confirmed by comp arison with the well-established and accurate results for about 500 asteroids provided by 
j 20121) and the n applied to determin ing H and G for the Pan-STARRS 1 asteroids using both the 


Pravec et al 


Muinonen et al. 


(120101) and 


Bowell et al 


(119891) phase functions. 
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1. Introduction 

Asteroid diameters are critical to understanding their dynamical and morphological evolution, potential 
as spacecraft targets, impact threat, and much more, yet most asteroid diameters are uncertain by > 50% 
because of the difficulties involved in calculating diameter from apparent brightness. The problem is that an 
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asteroid’s apparent brightness is a complicated function of the observing geometry, their irregular shapes, 
rotation phase, albedo, lack of atmosphere, and their rough, regolith-covered surfaces. Most of these data 
are unknown for most asteroids. The issue has been further confused because catalogued apparent magni¬ 
tudes for individual asteroids may have been reported by numerous observers and observatories over many 
years (even decades) in a variety of photometric systems with varying concern for ensuring accuracy and 
precision. This work describes our process for calculating asteroid absolute magnitudes (from which diam¬ 
eter is calculated) and their statistical and systematic uncertainties for hundreds of thousands of asteroids 
using sparse but accurate and precise data from a single observatory, the Pan-STARRS 1 facility on Maui, 
HI, USA. Our technique is suited to estimating absolute magnitude s when the phase curve coverage is even 


more sparse than those obtained by the Palomar Transient Factory (ILaw et al 


20091) . 


An asteroid’s absolute magnitude, H, is the apparent Johnson V-band magnitude, m, it would have if 
observed from the Sun at a distance of 1 au (i.e. observed at zero phase angle and 1 au distance). Accurate 
measurements of // as a function of time, together with infrared, polarimetric and radiometric observations, 
can provide crucial information about an asteroid’s size and shape, geometric albedo, surface properties and 
spin characteristics. 

In 1 985 the Internationa l Astronomical Union (lAU) adopted the two-parameter phase function devel¬ 
oped byO 


Bowell et al 


(119891 hereafter B89), OB(a; Hb, Gb), describing the behavior of the apparent magni¬ 


tude: 


m(r. A; Hb, Gb) = 5 log(rA) -i- 02j(a; Hb, Gb) 


( 1 ) 


where A represents the topocentric distance, r the heliocentric distance, and a{r. A) is the phase angle, the 
angle between the Earth and Sun as observed from the asteroid. We denote absolute magnitude in the B89 
system as Hb with a corresponding slope parame ter, Gb, that depends in a non-an alytical manner on (at 
least) an asteroid’s albedo and spectral type (B89: lLagerkvist and MagnussonLll99(lh . The slope parameter 
determines how strongly the apparent brightness of an asteroid depends on the phase angle and accounts for 
the properties of scattered light on the asteroid’s surfaces. Gb has an average value of ~ 0.15 (B89) for the 
most numerous S and C-class main belt asteroid taxonomies. An accurate determination of both Hb and Gb 
requires a wide and dense time coverage of the object’s apparent magnitude. Therefore, it is not surprising 
that only a few tens of slope parameters were measured before the advent of dedicated CCD asteroid surveys. 

The B89 phase function was very successful, but observations in the past twenty years have shown it can 
not reproduce the opposition brightening of E-type asteroids, the linear phase curve of the E-type asteroids. 
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and fa ils to accurately predic t the apparent brightness of asteroids at small phase angles. To address these 


Muinonen et al 


(1201(1 hereafter MIO) introduced an alternative phase function, (pM, with two slope 


parameters, Gi and G 2 that uses cubic splines to more accurately describe the behavior of the apparent 
magnitude. An alternative MIO formulation with a single slope parameter, G 12 that is denoted in our work 
as Gm, can be used when the data are not sufficient to derive the values of the two-parameter formulation i.e. 
m - 5 log(rA) -I- Their phase function was constructed such that ~ and the average 

asteroid would have a slope parameter of G^ ~ 0.5. This form of the phase function can provide better 
apparent magnitude predictions but der ivation of Hm and Gm st ill requires extensive light curve coverage 


and well-calibrated observational data (lOszkiewicz et al 


2 OI 2 I) . The IAU adopted the MIO (//,G1,G2) 


system as the new photometric system for asteroids in 2012. 

In the remainder of this work we use H and G to represent ‘generic’ absolute magnitudes and slope 
parameters respectively, and use the subscripts B and M on each parameter when referring to the values 
calculated using the B89 and MIO phase functions respectively. We implemented both functions to facilitate 
comparison with 1) past work that used the B89 parameterization and 2) future work that will use the now- 
standard MIO implementation. When we use Gm we specihcally mean the MIO G 12 parameter. 

The accuracy of most repo rted absolu te ma gnitud es is poor due to the lack of good photometry and 


limited phase curve coverage. 


.Turic et al 


(e.g. 


2002h first reported a systematic error of about 0.4 mags 


in the MFC’s absolute magnitudes which the MFC (and others) now attempt to address with observatory- 
dependent corrections to the reported apparent magnitudes. 

The determination of G has traditionally been even more of a challenge — they are so difficult to measure 
that they have only been calculated for «: 0.1% of asteroids and, even then, the uncertainty is usually large 


dFravec et al 


20121). An accurate measurement requires dense coverage of the phase curve and observations 


at different viewing aspects on the asteroid i.e. sub-solar positions. The vast majority of asteroids have no 
measured slope parameter so the average values of Gb - 0.15 or Gm - 0.5 are used. This assumption trans¬ 
lates into a systematic error in an individual asteroid’s H and G, and large uncertainty on the distribution of 
the parameters in the population. The problem is particularly acute for objects that have been observed only 


at large phase angles e.g. re sonant objects like 3753 Cruithne dde la Fuente Marcos and de la Fuente Marcos 


2013 


Wiegert et al 


1997h . and objects that orbit the Sun entirely within Earth’s orbit (IZavodnv et al 


2008) 


for which absolute magnitudes might be in error by up to about 1 mags. 

In summary, the problems with our current knowledge of asteroid absolute magnitudes and slope param¬ 
eters are due to; 
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1. Reporting observations to the Minor Planet Center (MPC) in non-standard filters and/or without accu¬ 
rate calibration. 

2. Not performing the color transformation from the hlter used for an observation to the Johnson V-band 
for an asteroid’s (usually unknown) color. 

3. The lack of information about the photometric uncertainty on each observation reported to the MPC 
so that it must be statistically ‘back-calculated’ for each observatory (or observer) from historical 
observations. 

4. The MPC database storing photometric values with only 0.1 mags precision. 

5. Assuming that Gb = 0.15 for all asteroids that do not have a reported value for the slope parameter. 

6. The accepted ‘standard’ average slope para meter of Gb = 0.15 for S and C class asteroids being 


different from the actual value of Gb - 0.20 (IPravec et al 


2012h . 


7. Sparse observations (in time). The lack of information about their rotation amplitudes induces an error 
and uncertainty in H. 


8. Selection effects (IJedicke et al 


2002h that bias the discovery of asteroids towards their rotation am¬ 


plitude maxima which induce a systematic error in their derived H. 

9. Most of the effort in deriving N and G focuses on their statistical uncertainties when the systematic 
uncertainties dominate. 

In this work we address each of these issues and derive the (Hb,Gb) and {Hm,Gm) parameters for 
known asteroids in the inner solar system out to, and includ ing, Jupiter’s Troian asteroids. All the data 


were acquired by a single wide-held survey, Pan-STARRSl (IKaiser et al 


20101) . in standard hlters with 


measured transformations to an accepted photometric system yielding photometric uncertainties that are 
typically about an order of magnitude smaller than earlier surveys. We use a Monte Carlo technique to 
measure the systematic errors introduced by hlter transformations for unknown spectral types, unknown G, 
and the unknown asteroid spin and amplitude. 

2. Pan-STARRSl asteroids. 


The Panoramic Survey Telescope and Rapid Response System’s prototype telescope (Pan-STARRSl; 


Kaiser et al. 


20101) was operated by the PS 1 Science Consortium during the time period in which the data 


used in this study was acquired. The telescope has a 1.4 gigapixel camera (ITonrv and Onaka . 


20091) and 


1.8 meter f/4 Ritchey-Chretien optical assembly and has been surveying the sky since the second half of 
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2011. Although the scientihc scope of the survey is wide — including the solar system, exoplanets, brown 
dwarfs, stellar astronomy, galaxies, cosmology, etc. — most of the data products are suitable for asteroid 
science. About 5% of the survey time was dedicated to the ‘Solar System’ (SS) survey (more accurately 
a survey for near-Earth objects, NEO) through the end of 2012, was increased to about 11% from then till 
2014 March 31, and the system is now 100% dedicated to NEO surveying. 

Pan-STARRS 1 surveys in six broadband bite rs, four of which were designed to be similar to the Sloan 
Digital Sky Survey photometric system (SDSS; 


Eukugita et al 


199611 . Most of the observing time was 


devoted to the Stt survey of the sky north of -30arcdeg declination for which each held was observed up 
to 20x/year in each of 5 biters — gpi, rpi, ipi, zpi and ypi. In the Stt survey the same held is observed 2 
or 4 times on a single night in 30-40 second exposures obtained within about an hour. The dedicated solar 
system survey used only the wide-hand wpi biter that is roughly equivalent to gpi + rpi + ipi with 45 second 
exposures and a cadence of ~ 20 min to image the same held 4x/night. The SS survey typically included 
belds within about 30 arcdeg of opposition or at small solar elongations ranging from 60 arcdeg to 90 arcdeg 
of the Sun. 

I mage processing was performed automatically and almost in real time by the I mage Processin g Pipeline 


(IPP; 


Magniei . 


2006h . Transient objects were identibed after ‘difference imaging’ (ILupton , 


2007h in which 


two consecutive images were convolved and suhbacted to identify moving, or stationary but variable, tar¬ 
gets. The photometric calibration until May 2012 was based on combined buxes of bright stars from Tycho, 
USNO-B and 2MASS catalogs. Since that time the entire northern sky has been imaged by Pan-STARRS 1 in 
all 5 biters allowing the development and use of the Pan-STARRS 1 st ar catalog with ‘ubercalibrated’ mag 


nitud es and zero points providing photometric uncertainties of ~ 1% (ISchlabv et al 


2012 ; 


20131). 


Magnier et al 


Moving tran sient detections are id entibed and linked into backlets by the Moving Object Processing 
System (MOP S; Denneau et al. 


(IMilani et al 


2013h and tracklets are associated with known asteroids by known server 


2008h . As of May 2015 the Pan-STARRS 1 MOPS team has submitted ~ 16,700,000 positions 


and magnitudes of 575,000 known asteroids to the MPC representing 85% of all numbered asteroids. During 
the same time period the system discovered ~41,000 asteroids, among them about 850 NEOs and 46 comets, 
and reported about 2,500,000 detections of unknown asteroids to the MPC. About ~ 42% of the detections 
were in the wpi biter acquired during the solar system survey while only about 9% were in the ypi and zpi 
bands. 

To ensure a consistent data set of high quality photometry (Pig. [U we restricted the detections used in 
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this study to known asteroids in the inner solar system (out to and including Jupiter’s Trojans) with multi¬ 
opposition orbits acquired during a sub-set of the Jin and solar system surveys between February 2011 and 
May 2012 (see Table T) The detections were selected from the IPP’s calibrated chip-stage PSF-fit photometry 


(ISchlaflv et al 


20121) and were required to be unsaturated, with sn >5, and not blended with stars or image 


artifacts. The Pan-STARRSl IPP never implemented the capability of fitting trailed asteroid detections, so 
we restricted our data sample to asteroids that trailed by less than 5 pixels during the exposure, equivalent 
to the typical PSF-width of ~ 1 arcsec. This limited the maximum rate of motion of the asteroids to about 
0.75 deg/day, excluding most NEOs and even fast-moving asteroids like Hungarias and Phocaeas on the 
inner edge of the main belt. Our strict criteria resulted in a set of more than one million detections of 
approximately 240,000 asteroids 


Table 1: Percentage of Pan-STARRSl asteroid detections in each filter in the time period from February 2011 to May 2012 (values do 
not add to 100% due to rounding). 


Band 

^pi 

'•pi 

/pi 

ypi 

^P1 

Wpi 

Fraction (%) 

18 

20 

17 

2.2 

6.2 

36 


Despite the enormous number of asteroid detections there are only about 10 detections/asteroid and 
each object is observed on average on only ~ 3 different nights over a phase angle range spanning about 
7 arcdeg (Fig. [TJ. Therefore, the survey pattern does not typically allow the determination of an asteroid’s 
spectral type, rotation amplitude or period. The detections have a mean+RMS photometric uncertainty of 
0.04 + 0.02mags and average+RMS visual magnitude of 19.8 +1.2 mags. The photometric uncertainty mode 
is ~ 0.02 mags corresponding to sm ~ 50 detections. This surprisingly high value is due to our selection 
criteria: the multi-opposition objects were identified in earlier surveys with smaller telescopes so they are 
typically brighter when observed with Pan-STARRSl. Note that only ~ 1% of the detections in our data 
sub-set have a photometric uncertainty greater than the 0.1 mags precision provided by the MPC. 

3. Method 

This work introduces a Monte Carlo technique to determine H (and G when possible) and its statisti- 
cal+systematic uncertainty based on the generation of synthetic asteroids (clones) that are each consistent 
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Photometric uncertainty (mag) Magnitude 


Figure 1: Characteristics of PS 1 asteroid detections used in this work, (clockwise from top left) number of detections per object, phase 
angle range per object, apparent V-magnitudes, and photometric uncertainties per detection. 

with the known asteroid. The clones explore the phase space of light curve rotation amplitudes, periods, 
colors and slope parameter in an attempt to replicate the observed apparent magnitudes. Each clone’s obser¬ 
vations are evaluated individually in the fitting process to derive H and G in the same manner as the actual 
observations so that the distribution of values for each object’s clones provide a measure of the systematic 
errors in the values. 
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3.1. Step 1: Initial fit for H and G 


The first step is essentially identical to the typical technique for calculating H and G; we fit the apparent 


the Levenberg-Marquardt least-squares fitting technique (ILevenberg 


1944 


Marauardl 


V-band magnitude to the B89 and MIO phase functions us ing the IDL procedure mpfit2dfur|j that employs 


1963) to minimize 


the variance between the detections’ apparent magnitudes and the values predicted by the models. We con¬ 
verted the Pan-STARRS 1 apparent magnitu des to V-band using taxo nomy-dependent filter transformations 


if the asteroid’s taxonomy was specified by 
(see Table|2|). 


Hasselmann et al 


(120121) and, if not, the mean S-)-C class color 


Table 2: Asteroid magnitude transformations from Pan-STARRS 1 AB filter magnitudes to the Johnson-Cousin V- system based on 


Tonrv et al 


' 2Q12h . Solar colors are also included for reference. 


Taxonomy 

V-gpi 

V-rpi 

V-ipi 

V-Zpi 

V-ypi 

V-Wpi 

Sun 

-0.217 

0.183 

0.293 

0.311 

0.311 

0.114 

Q 

-0.312 

0.252 

0.379 

0.238 

0.158 

0.156 

S 

-0.325 

0.275 

0.470 

0.416 

0.411 

0.199 

c 

-0.238 

0.194 

0.308 

0.320 

0.316 

0.120 

D 

-0.281 

0.246 

0.460 

0.551 

0.627 

0.191 

X 

-0.247 

0.207 

0.367 

0.419 

0.450 

0.146 

Mean (S-i-C) 

-0.28 

0.23 

0.39 

0.37 

0.36 

0.16 


The initial fits also use the mean class-depen dent G provided in Table|3]if the taxonomic class is specified 


in the SDSS database (IHasselmann et al 


2012 ) but, if the class is not known, we use the mean of the S - 


and C-class values: Gb - 0.15 (B89) and Gm - 0.53 (lOszkiewicz et al 


2012h respectively. 


The initial fits provided the absolute magnitudes in both photometric systems, Hs.i and Hm.i, that were 
the inputs to the next step in the pipeline. 


3.2. Step 2: Generating asteroid clones 

Our final H and G estimates are the result of Monte Carlo (MC) simulations that require the generation 
of synthetic ‘clones’ for each of the asteroids in our sample. Each of the clones is generated with its own 


’ Martiwardt IDL library, http://www.physics.wisc.edu/-craigm/idl 
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Table 3: Average slope parameters, Gb and Gm = 
Gi 2 , adopted in this wor k for 5 asteroid t axo¬ 
nomic classes as measured by[l 


Pravec et al 


'2012)and 


Oszkiewicz et al. 


' 20121) respectively. The 6th row 


provides ‘standard’ averages for the dominant S and 
C taxonomies. 


Taxonomic 

Class 

G = Gb 

+(RMS) 

Gi2 = Gm 

+(RMS) 

Q 

0.25+0.13 

0.41+0.14 

s 

0.24 ± 0.06 

0.41+0.16 

c 

0.15+0.09 

0.64 + 0.16 

D 

0.09 + 0.09 

0.47 + 0.14 

X 

0.20 + 0.09 

0.48 + 0.19 

S-hC 

0.15 

0.53 


color, slope parameter, rotation period, light curve amplitude and phase, where each of the parameters is 
selected from a unbiased distribution as described below. 


3.2.1. Clone colors 

Our pipeline can assign each clone the color of its parent asteroid (if known) or, when the parent’s color 
is not known, a random color based on an appropriate mix of taxono mies as a function of semi -major axis. 
About 16% of the asteroids in our sample have taxonomies dehned by|j 


Hasselmann et al 


(1201 2h (SDSS). 


We implemented this technique by dividing the inner solar system into 4 zones (see table |4]i; NEO-like 
(a < 2au), main belt {2au < a < 3.7 au), Hildas (3.7 au < a < 4.5 au) and Trojans (4.5 au < a < 6.0au). 
The semi-major limits dehning the zone edges were set at or near a minimum in the number distribution as 
a function of semi-major axis and by the availability of published taxonomic distributions. The exact values 
make little difference to this work. We used the publi shed, debiased taxonomic distributions in Table|4]in the 


4 zones with the qualihcation that for the main belt (iMothe-Diniz et al 


20031) we aggregated many related 


taxonomic types into 3 broad spectral classes: S-class=(A, AQ, AV, O, OV, S, SA, SO, SQ, SV, V, L, LA, 
LQ, LS), X-class=(X, XD, XL, XS), and Q-class=(Q, QO, QV). We required that the fraction, f(c,z), of 
asteroids with spectral class c in zone z satishes Yjcfic,z} - 1. I n the main belt, zone 2, w e were able to 


generate the taxonomies as a hner function of a as provided by (IMothe-Diniz et al 


2003h with a similar 


9 





















requirement that 2c /(c, a) = 1 at each semi-major axis. 


Table 4: Taxonomic distribution of asteroids in 4 semi-major axes zones used in this work. The main 
belt values are given below at a representative a = 2.5 au but we generated the clo ne taxonomies as a 
smooth function of semi-major axis in the range 2.0 au < a < 3.2 au as specified by 


Mothe-Diniz et al. 


200.3h. 


Taxonomy 

zone 1 

NEO-like" 

a <2au 

zone 2 

MB'’ 

a ~ 2.5 au 

zone 3 

Hilda" 

3.7 au < a < 4.5 au 

zone 4 

Trojans'' 

4.5 au < a < 6au 

Q 

14 

0 

0 

0 

s 

23 

61 

0 

0 

c 

10 

30 

7 

10 

D 

18 

0 

67 

80 

X 

35 

9 

26 

10 


Stuart and Binzel 

2004 


Mothe-Diniz et al. 

(200: 

i) 

Grav et al. 

(2012i 


Grav et al. 

(2012li 


3.2.2. Clone slope parameters 

We assigned slope parameters to the clones as a function of their assigned taxonomic class (c). i.e. 
the k‘^ clone was assigned a slope parameter Gk(c) = ran[G(c), crG(c)] where ran[x, y] is a random number 
generated from a normal distribution with mean x and standard deviation y, and G(c) and crc(T) are the 
mean and RMS of the distribution of slope parameters for class c, respectively (Table[3). 


3.2.3. Clone rotation periods, amplitudes and phases 

The sparse Pan-STARRSl data did not allow us to measure any asteroid’s rotation period and light 
curve amplitude. Furthermore, < 2% of the asteroids in our sample have measured ligh t curves reported 


in the asteroid light curve database (LCD 




Warner et al. 


2009; 


Waszczak et al 


2015 ) The lack of this 


^ The asteroid lightcurve database is publicly available at http://www.minorplanet. in£o/lightcurvedatabase.html 
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information introduces systematic uncertainty and error into the absolute magnitude and slope parameter 
determination. We quantified these effects using our Monte Carlo technique with synthetic sinusoidal light 
curves for each clone. 

Asteroid brightness variations on the hours-to-days timescales are usually caused by their non-spherical 
shape and rotation (the exceptions are for the unusual cases where the phase angle changes rapidly for 
close approaching NEOs, for multiple-systems in which brightness changes can occur if the objects transit 
or eclipse each other, and for objects with significant color variations). We assumed that the observing 
geometry (i.e. phase angle) effect on the asteroids’ light curves are negligible in the Pan-STARRSl data 
because of the limited range in phase angle coverage in our sample (Fig.[T]i. For the purpose of generating 
the clones’ light curves we assumed that all the objects are triaxial ellipsoids that generate simple sinusoidal 
light curves with peak-to-peak amplitude A, period P, and rotation phase 6. The offset from the unmodulated 
light curve at time t is then Am{t) - A sinilntlP + 0)12. 

Fight curve amplitudes tend to be larger for smaller asteroids (see Fig.|2] IWarner et al.ll2009h . probably 


because the smaller objects tend to be more irregularly shaped. Overall, the set of measured amplitudes and 
periods will be larger and shorter respectively than the true distribution because of observational selection 
effects, larger amplitudes and shorter periods are easier to detect and measure. 

To reduce the impact of t he light curve amplit ude and period selection effects we employed the debi- 


ased distributions derived by 


Masiero et al 


(l2009h that are representative of asteroids with ~ 18 (the 


average+RMS absolute magnitude in their study was 17.7 + 1.4mags). i.e. for objects with ~ 18 they 
provide the cumulative fraction of asteroids, Famp.MasiA), with light curve amplitudes < A. We empirically 
estimate the cumulative distribution of light curve amplitudes at other absolute magnitudes Famp{A,H) by 
‘normalizing’ to the median at 7/ = 18 from the median at other values; 

-'^me£/( 18 ) 


FampiA, H) = F, 


amp. Mas 


A X 


Affiedifi) 


( 2 ) 


where A^pAFl 

) is an 

Warner et al. 

2009) 


20091). Thus, given a clone’s initial ( ^3.Il l absolute magnitude, 77,, we generated a random 


light curve amplitude for the clone according to the cumulative fractional distribution given by eq. |2] 

We followed a similar procedure in assigning each clone a rotation rate R or, equivalently, a rotation 


period P = l/R. 


Masiero et al 


(120091) also provide the data from which we derive the cumulative fraction 


of asteroids, Frot,Mas(R), with rotation rates < R. Once again, their results are representative of asteroids 
with 77 ~ 18, about 2 mags fainter than the mean value in our data sample, so we developed an empirical 
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Figure 2: (left) Asteroid light curve amplitudes vs. absolute magnitude {Hb) from the LCDB IWamer et al 


20091) . The solid gray 


cui've represents the size-dependent moving median in 1.0 mags wide bins, (righ t) Measured asteroid spin rates (periods are provided 
on the right) vs. absolute mag nitude {Hb ) from the LCDB iWamer et al.L 120091) . The solid gray curve represents the size-dependent 
upper strength limit derived bv lHolsappld 120070 . 


technique to extend their cumulative fractional rotation rate distribution to other absolute magnitudes. 

Asteroids with diameters >100 meter (H < 23) have an empirically observed upper limit to their rotation 
rate of about 12rev/d av (Fig. |2]i and about 99% of the distribution of debiased spin rates are < 12rev/day 


(IMasiero et al 


2009j). Asteroids larger than a few tens of kilometers (H < 12) have an even more restricted 


upper limit to t heir rotation rates. W e empirically dehned an RmaxiH) as illustrated in hg.|2]and ‘compress’ 


or ‘expand’ the 
at any H: 


Masiero et al. 


(l2009h distribution as necessary to create the cumulative fractional distribution 


FrotiR, H) - Frot,A 


R X 




RmaxiH) 


(3) 


Once again, given a clone’s initial ( ^3.lb absolute magnitude, //,, we generated a random rotation rate for 
the clone according to the cumulative fractional distribution given by eq. [3 
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Finally, the rotational phase Ok for the A:'* clone was generated from a random uniform distribution in the 
range [0 arcdeg,360arcdeg). 

Our light curves were simple sinusoids even though we understand that real asteroid light curves can be 
much more complicated. The technique could easily be extended to incorporate actual light curve properties 
or a more realistic distribution but i) only a tiny fraction of known asteroids have measured light curves ii) we 
will show below that our results are not particularly sensitive to the actual light curve parameter distribution 
and iii) if the actual light curve is known then there is no need for any of the methods developed here. i.e. 
this method only applies to the 98% of asteroids that do not have measured light curves. Since this is a 
preliminary work we have not made any effort to remove those asteroids that have published light curves. 


3.3. Step 3: Refining H and G (First Monte Carlo simulation). 


The hrst Monte Carlo (MC) simulation yields our MC estimate for FI and G from the sparse Pan- 
STARRS1 phase curve coverage data. As described in detail above, we created 500 clones of each object 
where the k'* clone was assigned a taxonomic class (color) Ck, light curve amplitude A*, and period Pk. 
We then ht for each clone’s absolute magnitude, slope parameter and light curve phase, (//', G^, 6'^), by 
minimizing the;^^^ with respect to the actual observations: 

n2 

(4) 


A?, 


:,obs 'y j 

j=l L 


Mkitj; //', G', 6»p - m(tj) 


6m(tj) 


where n is the number of observations (detections) of the object, m(tj) is the actual object’s observed appar¬ 
ent magnitude, 6m(tj) is the reported uncertainty on the actual Pan-STARRS 1 apparent magnitude for that 
observation in the original hlter, and irik is the clone’s predicted apparent magnitude at the actual time of 
observation, tj, in the Pan-STARRS 1 filter in which the observation was made, with the clone’s appropriate 
color transformation (Amkitj)', Table|2]i: 


mk(tj) = 5 \og[r{tj)A(tj)] + <h[a(f)); H[, G^] H- Ak sin[2ntj/Pk + 0^]/2 H- Amk(tj), 


(5) 


and O is the B89 or MIO phase function. 

The ‘best’ clone is the one (k*) that produces the minimum and we adopt that clone’s and G^, 
values as our MC estimate for the object’s absolute magnitude and slope parameter. The process was run 
separately for both the B89 and MIO phase functions to provide ourMC estimates for (Hg, Gg) and (Hm, Gm) 
respectively. To avoid unphysical values the htting process required that -0.25 <Gg< 0.8 and -0.5 < Gm < 
1.5. 
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We found that 500 clones provides a good balance between the computation time and our ability to 
estimate the uncertainty on the absolute magnitudes and slope parameters. It is likely that when there are 
only a small number of detections that the number of clones could be decreased but we did not pursue this 


simplification. Whe n the number of det ections becomes very large then our technique becomes unnecessary 


as either traditional (IPravec et al 


20121) or sparse light curve fitting (IMuinonen et al 


201C; 


Law et al. 


20091) 


becomes more effective. 


3.4. Step 4: Estimating uncertainties and error on H and G (second Monte Carlo fit). 


We estimated the uncertainties and errors on and G^, by htting for the absolute magnitude and slope 
parameter with purely synthetic light curves generated from the clone with the best ht. i.e. we re-applied 
the same method as described in Step 3 03.3b except that we ht the clones to the best synthetic object rather 
than the real object (we continue to use the sub-script k to refer to clones but the clones used here are distinct 
from the clones used in the last step): 


2 

k.syn 


z 


i=\ *- 


Mkitj-, HI, G[, ep - 


-.2 


6mk>itj) 


(6) 


where dmkdtj) - 5m{tj), i.e. the uncertainty on the synthetic observation at time tj was set to the uncertainty 
on the actual observation at time tj. 

If we let X generically represent either H or G then the combined statistical-i-systematic uncertainty on 
X is the standard deviation of the clones’ X distribution: 




(7) 


where X' is the average value of X for all the synthetic objects’ clones. Similarly, the combined statisti- 
cal-Hsystematic error on X is the average error on the values for the synthetic clones: 



n 



( 8 ) 


3.5. Verification 

We verified our method with two independent sets of synthetic data generated from real Pan-STARRS 1 
data: 1) 10,000 randomly selected known Pan-STARRS 1 objects, most of them with sparse phase curve 
coverage and 2) the 1,000 known Pan-STARRS 1 objects with the best phase coverage. To have better control 
over assessing our method’s validity we generated photometric magnitudes and uncertainties with synthetic 
absolute magnitudes (Hg and Hm) and slope parameters (Gg and Gm) at each real time of observation with 
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the known object’s orbit. We then employed our pipeline to calculate each synthetic object’s H and G to 
measure the statistical and systematic errors induced by our technique. Moreover, we tested two different 
scenarios for a ssigni ng light cure amplitudes and periods to the clones: 1) the debiased dis tributions from 


Masiero et al. 


(120091) . 2) and the observed distributions from the LCDB (IWarner et al 


2009). 


The result is that for both synthetic populations (sparse and dense phase curve coverage) and for both 
light curve amplitude-period relations (debiased and observed) the difference between the generated syn¬ 
thetic values and the values returned by our method was normally distributed with zero mean. i.e. our 
technique correctly derives the H and G. Use of the debiased or observed amplitude and period distributions 
does not affect the derived H and G at the level of photometric accuracy and uncertainty of the Pan-STARRS 1 
data with its associated phase curve coverage i.e. does not cause any systematic errors. 


4. Results & Discussion 

4.1. Absolute magnitudes comparison with 
We think 


Pravec et al. 


mm 


Pravec et al 


(l2012l) ’s detailed light curve study of ~ 500 asteroids sets the standard in measur¬ 


ing asteroid photometric properties. They provided only Hb (it was before the adopti on of the new lAUstan- 
dard) but that value should be identical to Hm- Our results agree extremely well withy 


Pravec et al 


(120121) for 


the 347 objects that appear in both data sets (fig.|3]l. The mean differences of= -0.06+0.02 mags 
and Hm - Hs.pra - 0.02 ± 0.02 mags are consistent with zero to within 3cr and Icr respectively, with better 
agreement for the new lAU standard photometric system of MIO. The RMS of eac h distribution is 0.36 mags 
and 0.29 mags respectively, due to the quadratic combination of the errors in both IPravec et al.l (l2012h ’s and 
this work. 

The distribution of HB-Hg.pra is quasi-normally distributed (fig.[3]l with an RMS of 0.31 mags including 
a tail extending to Hp - Hp Pr„ < -1. Inte restingly, the difference between our initial fits with assumed slope 
parameter 1 ^3. lb andl] 


Pravec et al. 


(12012h . Hbj - HB.pra, is roughly normally distributed with a mean error of 


-0.06 ± 0.02 mags and RMS of 0.26 mags. Thus, the simple, traditional, fitting procedure with assumed G 
to our high-precision but sparse data produces comparable absolute magnitudes to the MC technique. The 
power of the MC technique lies in its ability to estimate the true statistical and systematic uncertainty in the 
absolute magnitude due to the unknown parameters in the analysis. 

Our absolute magnitudes calculated with the MIO phase function (Hm) are better behaved (fig.[3]l in the 
sense that the distribution is more normally distributed. The initial fit to the sparse data in the MIO system 
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Absolute Magnitude Uncertainty (5H) 



Pravec et al 


2012h using the B89 (left) 


Figure 3: (top) Absolute magnitudes from our study compared with 347 objects in common with 
and MIO (right) photometric systems. The dashed line shows the results of the traditional initial fit ( ^3.U and the solid line provides 
the results of the MC fit ( jj3.3t . (bottom-lef t) Uncertainties and (bottom-right) estimated systematic errors on absolute magnitudes from 
our study compared with those reported bylj 


Pravec et al. 


' 20121 ). 


provided absolute magnitudes with mean systematic errors of 0.00+0.02mags and cr ~ 0.26 mags compared 
to the MC technique with a mean error of 0.02 + 0.02 mags and cr ~ 0.28 mags. The good behavior of both 
the MC and initial fits with MIO that results in a normal error distribution leads us to the conclusion that it 
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is superior for the de termination of absol ute magnitudes even for sparse data samples. 


We also used the 


Pravec et al 


(120 12h values to test our technique ( ^3.41 i for establishing the uncertainty 


and error on our measured absolute magnitudes. Their technique allows excellent control of all the statistical 
and systematic uncertainties in the H calculation because they observed targets for more than a decade in 
a systematically controlled program and had 2 to 3 orders of magnitude more data per object. Thus, they 
report H uncertainties about 3x less than our uncertainties and we can compare our measured uncertainties 
{6H) to the RMS spread of H - Hpra, and our measured error estimates to its average (fig.[3]l. 

As stated earlier, the real power of the MC technique is its ability to estimate the statistical and systematic 
uncertainties on the derived H and G values. Our estimat ed absolute magnitude uncertainties (SHb', fig. |3j 
^3.4b for the asteroids that overlap the lPravec et al.l(l2012h data sample have the expected poissonian distri¬ 


bution with a mean of 6hb - 0.36 + 0.01 mags using the B89 phase function (fig.|3]l, comparable to the RMS 
of 0.37 + 0.02 mags for the distribution of the error in our measurement, Hb - HB^pra, as expected. Similarly, 
our mean estimated systematic error of AHb = 0.03 + 0.02 mags agrees with the actual systematic offset in 
the Hb - HB^pm distribution. We can compare our estimated uncertainties and systematic errors in the same 
manner for the MIO phase curve. Our estimated mean uncertainty, 5Hm - 0.26 + 0.01 mags, is consistent 
with RMS (Hm - Hs^pm) = 0.28 +0.02mags and our estimated systematic error of AHm = 0.00 ±0.02 mags. 


is consistent with (Hm - HB^pra) - 0.02 + 0.02 mags. 


The good agreement between our results and those of lPravec et al.l (12012h illustrates the utility of our MC 


technique at measuring an asteroid’s absolute magnitude and estimating the associated statistical+systematic 
uncertainty and any systematic bias, even for sparse data sets with limited phase angle coverage. Further¬ 
more, the nice behavior of our results with the MIO phase curve and the good agreement between our Hm 
and Hp^pra provides evidence that Hm ~ Hp when care is taken to ensure that the photometric data is of 
excellent quality. 


4.2. Absolute magnitudes 

Having established the utility of our technique on a well-controlled data set in the previous section we 
now employ it on all the asteroids in our selected Pan-STARRS 1 data sample. We were able to calculate the 
absolute magnitudes with combined statistical and systematic uncertainties for more than 240,000 asteroids 
spanning the range from 6.4 < H < 26.5 (fig. S. The ~ 20 mags range corresponds to about a factor of 
10, OOOx in the diameters of the objects and spans the inner solar system from the NEOs to Jupiter’s Trojan 
asteroids. Our sample represents ~ 38% of all known asteroids in that range as of February 2014, with the 
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highest completion of ~ 75% from 10.5 < H < 11.0. 




Absolute Magnitude Uncertainty (8H) 



Absolute Magnitude Error (AH) 


Figure 4: (left) Absolute magnitudes (Hm and Hg) of 248,457 asteroids, (center) Uncertainties and (right) estimated en'ors in the 
absolute magnitudes derived with our Monte Carlo method using the phase functions of (gray) B89 and (solid) MIO. 



Absolute Magnitude 


Figure 5: Absolute magnitude uncertainty as a function of H using the B89 (grey) and MIO (dark) methods. 


The mean uncertainties of 5Hb - 0.30 + 0.01 mags and 6Hm - 0.25 + 0.01 mags (Fig|4]l show that 
the new lAU photometric scheme of MIO is better than B89 for the sparse Pan-STARRSl data and phase 
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coverage but this conclusion mis-represents the full utility of the MIO technique. For one, the MIO system 
uncertainty is almost uniform with 6Hm ~ 0.24mags for the entire range \Q < H < 20 (hg.EJ. Second, 
even though the two techniques yield approximately the same uncertainties for the faintest objects for which 
the uncertainty is dominated by the measurement statistics (fig. H), the B89 method’s statistical uncertainty 
is ~ 0.35% larger for bright objects (10 < // < 14). 

The mean of our estimated statistical+systematic error using the MIO method, \AHm\ - 0.02+0.01 mags, 
is comparable to the B89 method, \IS.Hb\ - 0.01 ± 0.01 mags (hg.lD). The error in the absolute magnitude for 
each asteroid is less than the estimated uncertainty in ~ 62% of all the asteroids in our Hb sample and ~ 73% 
in our Hm sample. The RMS of the |A//b| and |A//mI errors respectively of ~ 0.35 mags and ~ 0.25 mags 
conhrms that the new lAU photometric system is an improvement over the earlier one and, furthermore, the 
shape of the error distribution is more reasonable for A//m than AHb (note the peak of A/Z^ is shifted by 
0.05 mags from zero but the A//m peak is near zero (fig. 0. 

Overall, there is almost no difference between our MIO and B89 ensemble results for Pan-STARRSl 
asteroids and the mean difference Hm - Hb is 0.03 ± 0.01 mags with RMS of 0.22 mags (fig.lHl. The mean 
difference between the initial ht solutions is Hmj - Hbj - 0.05 ± 0.01 mags. 




Absolute Magnitude Difference Absolute Magnitude Difference 

Figure 6: (left) Difference between the MIO and B89 absolute magnitudes for the MC and initial fit solutions, (right) Difference 
between MC and initial fit solutions for the absolute magnitude using the MIO and B89 methods. 

On the other hand, the utility of restricting H and G analyses to data derived from well-calibrated single- 
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survey data is easily illustrated by comparing the results of our technique to the MPC database values that 
do their best to incorporate data from multiple telescopes and observers over many decades. The MPC 
currently only publishes absolute magnitudes using the B89 phase function and there is a mean difference of 


Hm - Hb.mpc — 0.26 + 0.01 mags and Hb - Hb^mpc - 0.22 + 0.01 mags between our technique and the MPC 
values. The consistency between the mean differences is at least reassuring and the RMS spread in values 
is due to 1) the systematics introduced by the MPC’s procedure that incorporates apparent magnitudes from 
many different observatories in many different passbands and 2) the systematics introduced by our sparse 
light curv e coverage. Given th at we established in ^4.1l that our technique works well in comparison to the 
(120 12h values, our conclusion is that the error is due to the MPC’s incorporation of 


‘standard’ 


Pravec et al 


photometry from different sites a nd filters over a lo ng period of time. The error reported here is less than 
the ~ 0.4 mags value reported by|. 


Juricetal 


(l2002h . but since the time of that study the MPC database has 


been further populated by photometry from Pan-STARRS 1 and other large surveys with better photometric 
calibrations than previous surveys. Hence, it is unsurprising that the Hb,mpc values approach their correct 
values over time. 

Our calculated uncertainties are about 2x larger than reported byl( 


Oszkiewicz et al 


(120121) who employed 


the entire MPC catalog for their input photometry and provided Hm and Gm for 421,496 asteroids — al¬ 
most double our sample. For comparison with earlier works they also provided Hb and Gb- Their work 
was very difficult as it required calibrating and correcting the systematic problems intrinsic to the various 
observatories and observers that contributed the photometric data in multiple hlters, but offered the advan¬ 
tage of an extensive data set with wide time and phase angle coverage i.e. much like the MPC technique 


Oszkiewicz et al 


(l2012h 


described in the last paragraph. The systematic offset between our Hm values andy 
of Hm - Hjif osk = 0-3 3 + 0.01 mags (hg . |7l) is similar to the offset derived between our results and the MPC. 


Juric et al 


(l2002h and 


Pravec et al 


(120121) reported a systematic offset of about 0.38 mags to 0.5 mags 


between their c alculated absolut e magn itudes and the values reported by the MPC. Those values are in rough 


agreement with 


Waszczak et al 


(120 15h who reported Hb and Gb from over 54,000 asteroids observed in g 


and R-band with the Palomar Transient Factory (PTF). They measured a mean value of Rptf - Vmpc + 0.00 
which implies a systematic offset of ~ 0.4 mags in the MPC absolute magnitudes because the average V - R 
for asteroids is ~ 0.4. Our values (fig.© are consistent with the MPC for Hb <11 mags and Hb > 19 mags, 
i.e. within <0.1 mags of the MPC absolute magnitudes (their reported precision), but are systematically 
higher than the MPC absolute magnitudes for llmags<< Hb < 19 mags. i.e. our absolute magnitudes 
are systematically/a/nfer than reported by the MPC and this would translate directly into predicting fainter 
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Figure 7: The difference be t ween the Hg and Hm absolute magnitudes calculated in this work and (left) the Minor Planet Center 
and (right) lOszkiewicz et alj i2012n . We compare our Hm to the MPC Hg because the two photometric systems should yield similar 


absolute magnitudes (in theory). 


apparent magnitudes than the MPC and, similarly, suggesting that the objects are smaller than predicted by 
the use of the MPC absolute magnitudes. The systematic difference reaches a maximum of ~ 0.35 mags 
at Hb ~ 14 in agreement with the earlier studies. This magnitude offset has implications for developing 
observing programs, selecting objects for followup, and for studies of the asteroids’ size-frequency distribu¬ 
tion. 
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Figure 8: The difference between the absolute magnitude Hmpc reported by the MFC using the B89 phase function and this work’s 
Hb value as a function of absolute magnitude (< 100 random data points are shown at each discrete Hmpc value in order to reduce 
confusion). The thick solid black line represents the average in each 0.1 mags wide bin and the standai'd en'or on the mean is shown 
with en'or bai's. The error bars ai‘e about the width of the line for 13 < Hmpc <18. 
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4.3. Slope Parameters 

The vast majority of Pan-STARRS 1 asteroids offer only sparse phase angle coverage (Fig.lTJ for the de¬ 
termination of the slope parameter but our MC technique should provide a realistic estimate of the statistical 
uncertainty and systematic error when the phase angle coverage is not too large and the detections are not in 
multiple apparitions. 

The Gb distribution (fig. |9]l is very wide with a peak near 0.15, the default slope parameter for objects 
of unknown spectral class (most of the asteroids in our sample). The distribution is artificially constrained 
between the lower and upper limits (-0.25 <Gb < 0.8). On the other hand, the Gm distribution has a broad 
peak centered on Gm ~ 0.5 superimposed on a roughly flat distribution of slope parameters between our 
artificial limits (-0.5 < Gm < 1-5). The large peak near Gm - 0.2 that contains ~ 30% of all Gm values is 
due to a dis continuity in the MIO ph ase function, it is not an error in our implementation. In comparison. 


8% of the 


Oszkiewicz et al 


(120 12h Gm values were also ~ 0.2. Our technique is particularly sensitive to 


the function discontinuity and has a propensity to drive the fitted Gm value to 0.2 when the number of data 
points is small. We suggest that future attempts to use the MIO phase function flag and address this situation, 
perhaps by forcing Gm - 0.5 in those cases. 

The slope parameter uncertainty {6G) distributions have peaks at zero corresponding to the ~ 24% of 
cases in both methods where the MC technique did not converge and we fixed the slopes. Those Gb that 
were actually fit have a normal-like distribution with mean Gb - 0.18 + 0.01 and RMS of 0.05 (Fig. |9]l. 
Similarly, the Gm uncertainty has a normal-like distribution with mean at 0.29 + 0.01 and RMS of 0.17. The 
6Gm distribution is wider and shifted towards larger values than the 6Gb distribution because the Gm values 
are fundamentally larger than the corresponding Gb values. The percentage uncertainties ((5G/|G|, fig. 0 
in both slope parameters are very similar, suggesting that the two phase functions are equally effective for 
calculating the slope parameters, at least in the regime applicable to this data sample. The mean relative 
slope parameter uncertainties are ~ 34% and ~ 0.36% for Gb and Gm respectively, the large values being 
due mostly to the limited phase curve coverage. 

As expected, the slope parameter uncertainty depends on the phase angle coverage (Acr, fig. [Tol l. The un¬ 
certainty is artificially small at small phase angle ranges near zero because in these cases the slope parameter 
was mostly fixed at a pre-specified value. The uncertainty is largest for Aa ~ 5 arcdeg because at this phase 
angle range the slope parameter begins to be calculable, and the uncertainty drops at larger phase-angle 
ranges because the data provides stronger constraints on the shape of the phase function. However, even 
in the best case scenario, for phase angle ranges of > 30 arcdeg, the percentage uncertainty is still ~ 50% 
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Slope Parameter Error Relative uncertainty of G 


Figure 9: (top left) The slope parameters Gb (B89) and Gm (MIO) and (top right) their uncertainties and (bottom left) errors, (bottom 
right) Percentage uncertainty in the slope parameter {6G) with the (gray) B89 and (black) MIO photomeh'ic methods. 


for both phase functions. In any event, the number of objects in our data sample with large phase angle 
coverage is very small. Fig.[T0]also illustrates that the systematic errors introduced by our MC technique are 
not dependent on phase angle coverage. 


Pravec et al 


( 120121 ) provide acurate Gb slope parameters with uncertainties for more than 500 asteroids 


with densely covered light curves in a single pass band over a wide range of phase angles. The mean 


24 





























































Figure 10: (left) Average MC slope parameter uncertainty and (right) eiTor as a function of phase angle range using the (gray) B89 and 
(black) MIO phase functions. 


difference between this work’s Gb and Gb.pra is 0.00 + 0.02 with a ~ 0.28 for the 196 asteroids in common 
between the two data sets with derived slope parameters (Fig.[TT]i. The agreement between ourMC solution 


and the accurate work of 


Pravecetal 


( 120 12h using the B89 phase function suggests that our technique for 


calculating the slope parameter is viable for a large number of asteroids with sparsely sampled light curves. 


Furthermore, our technique allows us to estimate the mean error on the derived slope parameter, AGb = 
0.00 + 0.01, so the MC technique does not introduce a systematic bias. The mean statistical uncertain ty in 
the slope parameter for our data sample of 6Gb = 0.17 +0.01 is twice as large as thej] 


Pravec et al 


(1200 data 


set of 6Gb - 0.09 + 0.01 which could be interpreted as either surprisingly good, given the small number of 
observations and phase curve coverage of our data sample, or as an indication that measuring Gb is difficult 
even with a very good da ta sample. 


As described earlier. 


Oszkiewicz et al 


(120 12h derived asteroid slope parameters from photometry re¬ 


ported to the MFC from multiple observatories that used different filters and reference catalogs. They also 
had to deal with the fact that the MFC observation submission format did not allow reporting of photo¬ 
metric uncertainties. To reduce some of the associated problems they statistically calibrated the disparate 
datasets and used photometry only from major surveys. After excluding the artificial peak near Gm - 0.2 
(i.e., excluding the range 0.18 < Gm < 0.22), and including only those objects for which G was actually 
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Slope Parameter Difference (Gg-G^pp^) 



Figure 11: 
parameter 


uncertainty distributions for the same 196 objects for (solid) our MC values and (dashed) 


Pravec et al. 

20121 

Pravec et alj 

2012) 


Table 5: Mean slope parameters ± standard deviation {Gb, 
B89) derived in this work (PSl, second column) and by 


Pravec et al 


20121) (PRA12, third column) for the same ob¬ 


jects in 4 major taxonomic classes. The last column is the 
number of common objects that have a iHasselmann et al. 


' 2012h spectral classification (no D type asteroids satisfied 
our requirements on taxonomic identification and slope pa¬ 
rameter determination). 


Taxonomic 

Class 

Gb 

PSl 

Gb 

PRA12 

N 

Q 

0.11+0.16 

0.19 + 0.10 

3 

s 

0.16+0.26 

0.23 + 0.05 

32 

c 

0.03+0.10 

0.13 + 0.01 

4 

D 

n/a 

n/a 

0 

X 

0.21+0.30 

0.20 + 0.10 

9 
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Slope Parameter Difference 




Slope Parameter Difference 



Oszkiewicz et al. 


Oszkiewicz et al 


201 


i) 


Oszkiewic£_et al.l <2012h . (bottom right) Slope 


I2OI2I) . (bottom 


Figure 12: (top left) Difference between our MC Gb values (MIO) and 133,885 objects in com mon with Oszkiemcz_etalJ 42Q12h 
(top right) Slope parameter uncertainties for the same objects as determined in this work and by 
left) Difference between our MC Gm values (MIO) and 80,756 objects in com mon with 
parameter uncertainties for the same objects as determined in this work and by 


fit, there were 80,756 objects in common with our Gm values and 133,884 objects for comparison with our 
Gb- The wide and oddly-shaped d istribution of the difference in slope parameters between our MC tech¬ 


nique and 


Oszkiewicz et al 


(12012h (fig. [TZt illustrates the difficulty and large uncertainty in measuring G. 
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Table 6: Slope parameters derived in this work (PSl: Gb, second column; Gm fifth column) and by 


Oszkiewicz et al. 


' 2012h (OSK12: Gb, third column; Gm, sixth column) for the same objects in five differ¬ 


ent spectral classes. The forth and last columns are the number of objects in common between the two data sets 


with SDSS spectral classification iCarvano et al 


20101 . 


Taxonomic 

Class 

Gb 

PSl 

Gb 

OSK12 

N 

Gm 

PSl 

Gm 

OSK12 

N 

Q 

0.21+0.28 

0.20 + 0.10 

1324 

0.46+0.53 

0.54 + 0.22 

886 

s 

0.22+0.28 

0.19 + 0.22 

14686 

0.47+0.53 

0.55 + 0.20 

10231 

c 

0.18+0.28 

0.16 + 0.10 

7892 

0.58+0.55 

0.66 + 0.23 

5150 

D 

0.23+0.29 

0.19 + 0.12 

1321 

0.42+0.52 

0.61+0.25 

852 

X 

0.19+0.28 

0.18 + 0.11 

2073 

0.53+0.54 

0.59 + 0.24 

1428 


The distribution peaks at zero for the B89 phase function with (Gb - Gb,osk - 0.00 + 0.01) but there is a 
significant offset using the MIO phase function of Gm — Gm,osk - -0.06 + 0.01 (fig. US. The RMS of the 
difference is larger using the MIO (0.58) than with the B89 phase function (0.35) but this is expected due to 
the numerically larger expected values of Gm ~ 0.5. 

Fig.[l2]also illustr ates that our MC techniqu e yields slope parameters that are comparable or marginally 


better than the work of 


Oszkiewicz et al 


(120121) . even though our data sample includes much less photomet¬ 


ric data per object over a narrower phase angle range, presumably because of the Pan-STARRSl system’s 
superior photometry an d the use of measured ph otometric uncertainties. The mean uncertainty for 80,756 


objects in common with 


Oszkiewicz et al 


(12012h is 0.33 + 0.01 (RMS- 0.14 ) with our MC technique and is 


0.39 ± 0.01 (RMS= 0.18) for the values report ed bvlOszkiewicz et al 
_ Sl ope parameters are t axonomy-dependent (IHarris 


1989 


( 2012 ). 


2012 ; 


Pravec et al 


Lagerkvist and Magnusson , 


199C; 


Oszkiewicz et al. 


2012h but most of the objects in our Pan-STARRS 1 data sample are fa inter than known 


astero ids with well established taxonomies, so we relied on the SDSS spectral classification (iHasselmann et al 


with 


2012) to assess our meth od’s ability to dete ct the taxonomic-depend ence. We found 48 asteroids in common 


Pravec et al 


mm and 18,541 with 


Oszkiewicz et al 


(12012h (excluding values around Gm ~ 0.20) 


for which we cou l d com pare our calculated slope parameters. Our mean+RMS Gb values are consistent 
(12012h (Table l43T l but our RMS distribution is much larger and the common numbe r of 


with 


Pravec et al 


asteroids is very low. Similarly, our Gb and Gm values (Table |6]l are consistent 


Oszkiewicz et al 


(12012h but 


the RMS is distributions are large in both cases. There is a formal difference between the means of some of 
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the taxonomic classes but we do not consider them further because of the large uncertainties on each value 
and the large RMS of each taxonomic class’ G distribution. 

As discussed above in section 13.11 the phase curve coefficients Gb and Gm are functions of asteroid 
composition. Given the compositional trends of the inner belt being dominated by silicate S-type asteroids 
and carbon/volatile-rich asteroids in the outer belt, we sh ould expect to see these trends reflected in our 


derived phase functions. A similar study was performed by 


Oszkiewicz et al. 


(12012h in their analysis of the 


MFC database. They found correlations between their measured G12 and orbital elements throughout the 
main belt, reflecting the general compositional gradient and family structure. To study this in our database 
we selected the 51,864 asteroids with orbital semi-major axes 2.1 < a < 3.3 AU where the range of phase 
angles observed was Aa > 5° and there were N > 6 observations. We then calculated the running median 
values Gb and Gm as a function of orbital a over a range Aa - 0.05 AU. 




Figure 13: Moving average of Gm (top) and Gb (bottom) as a function semi-major axis. 


Figure [13] clearly shows clear a negative trend in Gb and a positive trend in Gm with orbital a. As 
Gb is larger for S-type than C-type asteroids, while Gm becomes smaller, this agrees with the established 
compositional gradient in the main-belt. For modelling purposes, these trends may be approximated by the 
relationships Gb = -0.103a -i- 0.446 and Gm = Q.237a - 0.175 within the main belt. The largest deviations 
from these relationships occur at the 3:1 Kirkwood gap at 2.50 AU, and at the 7:3 gap at 2.95 AU. This latter 
position marks where the S -type asteroids of the dom inant Koronis family of gives way to the T/X/K/D-type 


asteroids of the Eos family (IMothe-Diniz et al 


2005h. We note that the overall observed scatter in individual 


values is dominated by AG, although it will also be partly due to the large amount of compositional mixing 
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present in the main belt (IDeMeo and Carry 


20131). 


5. Availability 


The Pan-STARRS 1 absolute magnitudes and slope parameters with associated uncertainties as described 
herein are available on-line (Appendix Appendix A| ). The eventual goal is that the catalog will be updated 
with all the data from the entire 3 year Pan-STARRS 1 mission and then updated regularly with new data 
from the ongoing extended mission that is purely focused on the solar system. This effort will provide 
almost complete coverage of all known asteroids with extensive phase angle coverage and good number of 
detections per object. 


6. Conclusions 

Our work introduces a Monte Carlo method for calculating absolute magnitudes (H) and slope parame¬ 
ters (G) and their statistical uncertainties and systematic errors that is applicable to single apparition asteroid 
observations and designed to handle limited photometric data over a restricted phase angle ran ge. The tech 


nique’ s utility was confirmed by comparing our H and G values to the well-established results of lPravec et al 


(12012h for a limited number of objects. We then applied it to derive H and G with statistical uncertainties 


and systematic errors for ~ 240,000 numbered asteroids observed in the first 15 months of Pan-STARRS 1 ’s 
3-year nominal mission. The single-survey data, consistent image processing, and well-defined photometric 
calibration, eliminates many of the problems encountered in past attempts to measure absolute magnitudes 


We find that the 


Muinonen et al. 

(2010 

) phase function provides better results than the 

Bowell et al. 


(119891) phase function in terms of reducing the statistical uncertainty and systematic error on the absolute 


magnitude — both crucial to accurately predicting ephemeris apparent magnitudes and calculating asteroid 
albedos from H and measured asteroid diameters. There is a systematic //-dependent offset between the 
Minor Planet Center’s reported absolute magnitude and H derived in this work with a maximum offset of 
about 0.25 mags at // ~ 14. 

The measured slope parameters are generally in agreement with the results of [j 


Pravec et al 


Oszkiewicz et al. 


(120121) and 


(12012h but the statistical uncertainty and systematic error on any individual asteroid’s G is 


large due to poor temporal and phase-space coverage. 
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Appendix A. Pan-STARRSl asteroid database 


Version 1.0 of the Pan-STARRSl asteroid database is available at http: //www. ifa.hawaii . edu/NEO/ 
It provides derived H and G values for 248,457 asteroids with a total of 1,242,282 detections spanning the 
time interval from February 2011 to May 2012 as described in this work. The 18 column data hie is comma- 
delimited and each line represents a single asteroid. The columns are described in table lA.7l 
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Table A.7: Pan-STARRSl asteroid database vl.O column descriptions. 


Col. 

# 

Col. 

value 

Description 

1 

ID 

The object’s designation in the MFC's 5-character format. The MFC database is accessible online.“ 

2 

class 

The obiect’s taxonomic class as soecified bv the Sloan DipitalSkv Survev tHasselmann el al.. 20121 from the Flanetarv Data Sv.stem. version 1.1. available online^. 

NULL if unknown. 

3 

N 

number of detections used in the fit 

4 

Aa 

phase angle range 

5 

^B.i 

initial estimate of the absolute magnitude using the B89 phase curve 

6 

Hb 

absolute magnitude derived using our MC technique in the B89 photometric system 

7 

SHq 

uncertainty on the absolute magnitude in col. 6 

S 

AHb 

estimated error on the absolute magnitude in col. 6 

9 


initial estimate of the ab.solute magnitude using the MIO phase curve 

10 


ab.solule magnitude derived using our MC technique with the MIO pha.se curve 

11 

SHm 

uncertainty on the absolute magnitude in col. 10 

12 

AHm 

estimated error on the absolute magnitude in col. 10 

13 

Gb 

slope parameter derived using our MC technique in the B89 photometric system 

14 

8Gb 

uncertainty on the slope parameter in col. 13 

15 

AGq 

estimated error on the slope parameter in col. 13 

16 

Gm 

slope parameter derived using our MC technique in the MIO photometric system 

17 

^Gm 

uncertainty on the slope parameter in col. 16 

18 

AG^ 

estimated error on the slope parameter in col. 16 


“http://www.minorplanetcenter.net/iau/MPCORB/MPCORB.DAT 
*http://sbn.psi.edu/pds/resource/sdsstax.html 
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